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Abstract 

The heat capacity Cp of an array of independent ahiminum rings has been measured under 
an external magnetic field H using highly sensitive ac-calorimetry based on a silicon mem- 
brane sensor. Each superconducting vortex entrance induces a phase transition and a heat 
capacity jump and hence Cp oscillates with H. This oscillatory and non-stationary behaviour 
measured versus the magnetic field has been studied using the Wigner-Ville distribution (a 
time- frequency representation). It is found that the periodicity of the heat capacity oscilla- 
tions varies significantly with the magnetic field; the evolution of the period also depends on 
the sweeping direction of the field. This can be attributed to a different behavior between 
expulsion and penetration of vortices into the rings. A variation of more than 15% of the pe- 
riodicity of the heat capacity jumps is observed as the magnetic field is varied. A description 
of this phenomenon is given using an analytical solution of the Ginzburg-Landau equations 
of superconductivity. 
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1 Introduction 



Thanks to recent evolutions in micro and nanofabrication processes, new devices are emerging which 
permit the study of the specific heat of mesoscopic systems even for samples having very small mass (typ- 
ically under 100 nanogram) [U [5J |2l HI [3 [5] ■ This kind of innovative sensor opens up a great possibility of 
studies which can be carried out on nanosystcms. It enables the study of the phase diagram and the phase 
transitions specific to systems of small size. Indeed, the reduction of the dimensions of superconducting 
or magnetic systems leads to the appearance of new phase transitions J\ , or mechanisms proper to quan- 
tum physics which can be studied from a thermal point of view [SI El [To]. More specifically, the study of 
vortex matter in nano-engineered arrays of mesoscopic superconducting systems is of particular interest. 
The effect of the geometry, the thickness of the film or the size of the system compared to a physical 
characteristic length scale can be studied from a thermal point of view in order to better understand the 
physics of vortex in extreme or particular conditions. In general calorimetry has the advantage of being 
sensitive to any energy level in the system. Magnetization for instance is sensitive only to the magnetic 
aspect of the sample. Through calorimetric measurements any phase transition whatever its origin is 
expected to exhibit a signature. 

In a superconducting ring a supercurrent appears to screen or enhance a perpendicular external 
magnetic field (described by a vector potential A), so as the fiuxoid <j)' is fixed to an integer multiple of 
^0 — h/2e (the superconducting quantum flux); (p' is then given by: 



where m* , e* and iJl are respectively the mass, the charge and the velocity of the supercurrent carriers, 
(p is the magnetic flux threading the integration contour and n is an integer. The integer number n of 
fluxoid quanta in the ring is a good quantum number and can be seen as the number of vortices threading 
the ring, each vortex carrying a fluxoid quantum ^q. If this kind of systems has been widely studied 
theoretically [TTl [T^] and experimentally from an electrical point of view [TSl [Uj , by magnetic deco- 
ration [16] or through magnetization measurements [171 118j , little is known about the thermal behaviour 
of superconducting mesoscopic loops[31[51[Tni[ini[2I]- Here we report highly sensitive heat capacity mea- 
surements performed on an array of independent mesoscopic superconducting rings of size comparable to 
the superconducting coherence length ^(T) under an applied magnetic field. 

In a previous paper (see Ref. 0]), we have already demonstrated that multiple phase transitions be- 
tween states with different vorticities are accompagnied by discontinuities of the heat capacity. Each 
vortex entrance (or expulsion) is associated to a mesoscopic phase transition from the n to the n + I 
(or n — 1) state, and as the magnetic field is increased, an oscillating heat capacity is measured with 
a periodicity corresponding to a flux quantum 0o threading a loop. Under speciflc conditions (lowest 
temperatures, zero fleld cooling) it was also shown that several vortices could enter or exit the loop at 
the same time, leading to oscillations of the heat capacity with a periodicity of 200 at 0.85 K and 300 
at 0.70 K. An important feature distinguishing this article from Ref. |1] is that here we focus only on 
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the 100 periodic oscillations of the heat capacity: even at lowest temperatures the system is prepared 
to have only a 10o periodic component. Indeed as it was already nientionned in Ref.[4J, the magnetic 
history of the loops is very important and determines the way the sample behaves. If the magnetic field 
is swept from -40 mT down to zero even at 0.6 K the modulation of the heat capacity is (j)o periodic. 
The 300 periodic signal as shown in Ref.[4] is evidenced only when the loops are cooled in zero field 
before the magnetic field is swept from zero to a high value of H. This hysteresis of the heat capacity of 
mesoscopic sample has been observed on many different situations; this behaviour will be discussed in a 
forthcoming paper. Here we avoid the appearance of 200 and 30o periodic oscillations, in order to study 
the fine evolution of the periodicity of the heat capacity jumps with regard to the varying magnetic field 
and to its sweeping direction. The advantages of working at low temperature are the enhancement of the 
signal-to-noise ratio and the increase of the total number of vortices a loop can host. We stress here that 
in Ref. [4] we studied the temperature and magnetic field ranges enabling the appearance of multiquanta 
transitions. It was shown that the pseudoperiod of the heat capacity versus magnetic flux C{H) can take 
discrete values nx (j)o. In the present paper we focus on the 10o regime and show that the pseudoperiod 
of C{H) is continuously changing as H is swept. 

2 Experimental results and analysis 
2.1 Heat capacity measurements 

The sample studied in this work is composed of an array of A'^ = 4.5 x 10^ identical non interacting 
superconducting alumimun square loops (see the inset of Fig. [2j each loop has 2 /im side, w = 230 nm 
arm width, d = 40 nm thickness, for a total mass m = 80 ng of aluminum). The separation of neighboring 
loops is 2 /im. The inter-loop interaction is neglected; we indeed do not expect it to affect our results 
significantly because for our sample the mutual inductance (~ —20 fH) is much smaller than the self 
inductance (~ 5 pH). Thus the mutual magnetic flux always remains much smaller than 0o, and the 
energy of magnetic interaction between loops is much smaller than the free energy of a single loop. We 
will consider in the following that all the thermal signals are additive, hence the measured heat capacity 
is N times the heat capacity of a single loop. 

The N mesoscopic square loops are patterned by electron beam lithography on a home-made specific 
heat sensor and aluminum is deposited by thermal evaporation. The sensor [2j|4j is composed of a large (4 
mmx4 mm) and very thin (5 fim) silicon membrane suspended by twelve silicon arms. On this membrane, 
a copper heater and a NbN thermometer are deposited through regular photolithography. Resistances of 
these thin film transducers are measured by a four point probe technique. In the case of the Cu heater 
we can thus measure the power injected in the calorimeter. In the case of the thermometer the four wire 
measurement reduces the measuring noise, enabling the read out of very small temperature variations 
on which our specific heat measurement is based. Indeed using this device the total heat capacity can 
be measured by ac-calorimetry 22 . The principle of this technique lies in applying a sinusoidal current 
in the heater at the frequency /, resulting in oscillations of the temperature of the membrane at twice 
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this frequency. Thus, when measuring the voltage response on the dc-biased thermometer one gets the 
ampHtude (5Tac of the oscillations of the temperature and hence information about the thermal properties 
(thermal conductivity and heat capacity) of the silicon membrane and the nanosystems it contains. This 
method has been largely described in numerous publications [U [22l [H [2] . The calorimetry setup is cooled 
down to 0.55 K using a ^He cryostat. 

The correct working frequency is obtained by measuring the response function / x STg^c of the sensor. 
When this response is extremal, i.e. becomes independent of the frequency /, the system can be considered 
as quasiadiabatic and then i^Tac is only related to the heat capacity and to known parameters (frequency 
/ and injected power) [22] . We define the adiabatic plateau as the frequencies interval such as / x (5Tac is 
greater than 99% of its maximum. In our case we obtain / € [108; 149 Hz] at 0.6 K; in these conditions 
the thermal excitation is faster than the characteristic time of the heat loss to the thermal bath but slower 
than the heat diffusion time in the system (sensor and sample). The major advantage of working with 
lithographied loops is the high thermal contact between the nano-objects and the silicon sensor. In this 
quasiadiabatic limit, all superconducting loops can be considered as being at the same temperature. The 
amplitude (JTac of the oscillations of temperature can be tuned between 1 mK to 15 mK depending on 
the working temperature and on the resolution needed for the measurements. The ac-calorimetry enables 
averaging of the measured signal, thus reducing the error bar for each measured data point. Typically by 
averaging over 10 seconds, this apparatus allows measurements of heat capacity within 10 femto- Joule 
per Kelvin, which corresponds to an energy sensitivity as small as few atto-Joule (10^^^ Joule). 

The sensor is placed in the center of a large superconducting coil so that a tunable magnetic field H 
can be applied perpendicular to the plane of the loops. The area containing the loops is 2.5 mmx2.8 
mm large, and the difference between the field on the axis and the field at the border of the sample has 
been calculated to be less than 0.1%. This is regarded as negligible and hence the field is considered as 
homogenous all over the sample surface. A heat capacity measurement at 0.75 K is shown on Fig. [TJ 
The measurement is performed under a perpendicular magnetic field H swept from -45 mT to 45 mT at 
constant temperature. Two major signatures of second order phase transition from the normal phase to 
the superconducting phase can be seen at -35 mT and 35 mT through a heat capacity jump of 0.75 pJ/K. 
Apart from these significant changes in the heat capacity, subsignatures appear at intermediate magnetic 
fields. As it is expanded in the insets of the Fig. [l] oscillations of the heat capacity can be observed with 
a period of approximately 0.6 mT. This magnetic field interval corresponds to one superconducting fiux 
quantum 0o {(j)o ~ 5^) threading a square of 1.86 /im side, a contour which is included in the volume of a 
single aluminum loop (see the inset of Fig. [2]). As it was demonstrated earlier [4J, each vortex entrance (or 
expulsion) changes the vorticity n by one, leading to a phase transition characterized by a heat capacity 
jump. As the flux through a loop is swept such phase transitions occur with a periodicity close to (j)o and 
hence an oscillatory behavior appears. The amplitude of these oscillations (w 20 fJ/K) is of the order 
of the measuring noise (« 10 fJ/K by averaging over 10 seconds). So the oscillations are hidden in the 
noise when looking directly at the raw data of Fig. [1] Oscillations appear more clearly (as shown in the 
insets of Fig. [T|) by smoothing the raw signal through a 10-points adjacent-averaging process. 
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In order to better exhibit the osciUatory behavior, we perform a filtering of the background of the 
signal (general trend corresponding to low frequencies); then we calculate a Fourier transform of the 
oscillating signal isolated in this way. The result is presented in Fig. [2] where the modulus of the Fourier 
transform shows the oscillating contribution of the superconducting loops inside the superconducting 
area (H = -30 mT to H = 30 mT). A large peak is exhibited at the frequency v — 1.65 mT^^, but 
not much information can be extracted from this broad peak. Indeed since our aim is to study how 
vortices entrances or expulsions affect the period of the signal as the field is swept, this single spectral 
representation is not adapted. In fact, to extract relevant information of a non stationary signal, one 
needs to have a signal processing tool with a high resolution in both magnetic field H and frequency 
domain [v). Hence we will use a specific time- frequency representation, the Wigner-Ville distribution. 

2.2 Signal processing 

The magnetic field interval between two successive heat capacity jumps seems to evolve as the field 
is increased. This new feature, not mentioned in Ref. [4], is quite interesting if one wants to deeper 
understand the specific character of the metastability of the thermodynamic state in a superconducting 
ring. From electrical measurements on a mesoscopic superconducting disk, Baelus et al. have shown that 
the penetration and expulsion fields can depend on the measurement temperature 23 . Here we show that 
the penetration and expulsion fields depend on the magnetic history of the sample as well, i.e. they depend 
on the applied magnetic field and on its sweeping direction. In order to study how the periodicity of the 
Cp{H) curves evolves as the magnetic field is varied, we use the Wigner-Ville distribution (WVD)[24j, a 
time-frequency processing tool adapted to study non-stationary signals and giving a better resolution in 
frequency than a regular Fourier transform, especially when the characteristic frequency of the signal is 
changing continuously with time. For our data the role of time is played by the magnetic field. 

The main purpose of the WVD is to image the energy distribution of a non-stationary signal in the 
time-frequency space. Because any tool based on simple Fourier transform needs to find an appropriate 
balance between time and frequency resolution, the WVD solves this problem by calculating 

/■ + °° T T 

Wc,{H, j ^ C,{H -f -)C,{H - -)e-2— dr (2) 

where Cp{H) is the heat capacity measured versus magnetic field. The WVD gives the best compromise 
between temporal and frequency resolutions. However the WVD processing method has a disadvantage: 
the distribution is perturbed by the presence of cross-terms. These cross-terms emerge from interferences 
between characteristic frequencies. To resolve that problem, instead of using the raw definition given in 
Eq. (21) we use the smoothed pseudo Wigner-Ville distribution (SPWVD) given by 



5c,(i?,.)=mMi)r>< 

j+^g{s~H)C,{s + r)c;{s-i)ds 



(3) 



where h(T) is a short time window allowing a spectral smoothing to reduce cross-terms between shifted 
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time terms and g{s — H) is a frequency window allowing a time smoothing to reduce cross-terms between 
shifted frequency terms. Several kinds of windows can be used (see Ref. ^25j). In the same way as in 
classical spectral analysis, boxcar window function are prohibited, so we use in this work smoothing in 
time and frequency through Hanning windows |25j . By applying two Hanning windows (h and g) one 
gets the same result as a 2D convolution of the WVD by a 2D window. Consequently SPWVD causes 
a slight loss of resolution for the characteristic frequency but enables a strong attenuation of the cross 
terms. 

We use the SPWVD to analyse our Cp{H) signal. An example is shown in Fig. [3] This is a density 
plot of |5| in the H — v plane; this SPWVD is calculated from the Cp{H) data presented on Fig. [1] At a 
given magnetic field H (horizontal axis) the amount of power contained in the harmonic with frequency 
v (read on the vertical axis) is given by a colour code, from the deep blue for the lowest amplitudes 
to the red for the highest ones. No quantitative physics will be extracted from the amplitudes given 
by the color code, which is only used to localize the characteristic frequency read on the vertical axis. 
Roughly speaking, for a given field H, the ordinate of the red area is the local frequency of the signal 
in the vicinity of H, i.e. the inverse of the local periodicity of Cp{H). Thus we clearly get the evolution 
of the periodicity in the vortex expulsion area (negative magnetic field) and in the vortex penetration 
area (positive magnetic field): at -25 mT the local frequency is (1.50±0.04) mT^^ (i.e. the period is 0.67 
mT) and at 25 mT the frequency reaches (1.78±0.02) mT^^ (period — 0.56 mT). The change of period 
of the heat capacity oscillations is thus above 15%. This is an illustration that all the jumps in the heat 
capacity observed in that kind of experiment are not directly related to one (j)o, but are linked to the 
metastability of the thermodynamic states which will change with the applied magnetic field, and so with 
the number n of vortices in the loops. Either the periodicity of the jumps is not exactly an integer number 
of quantum magnetic flux (j)Q, or the contour through which the fluxoid is strictly quantized evolves with 
the magnetic field. This will be further discussed in the theory part of the article. 

The non-stationarity of the signal does not depend on the temperature; closer to Tc the same obser- 
vation can be made. It does not depend neither on the geometry because the same results have been 
obtained on circular superconducting rings of 1 fim and 2 fim in diameter. The crucial parameter is the 
sweeping direction of the applied magnetic field. If the absolute value \H\ of the field is decreased, the 
sample is in the situation where vortices are expelled outside the loops. The opposite behavior is obtained 
when vortices are penetrating the loops when \H\ is decreased. We clearly illustrate that point in Fig. |4] 
and Fig. O (the analogues of Fig. [T] and Fig. [3]) where the magnetic field is decreased from 40 mT. In 
this case, as the field is decreased, the frequency of the Cp jumps is increasing. This evidences the fact 
that the sweeping orientation of the field determines the evolution of that periodicity. The magnetic field 
sweep rate was the same («4.5 mT per hour) when getting the data points of Figs. [T]and[H and it can 
be varied between 1 mT/hour to 20 mT/hour without affecting the results. 
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3 Theory 



A reasonably good description of our experimental observations can be obtained in the framework of the 
Ginzburg-Landau (GL) theory of superconductivity [26j . In this theory the free-energy density in the 
superconducting state / is written as: 



/(f) = /„o™ai(r-) + a{T) mf)\' + ^ mft 



+ 



2m* 



-V A ijjir) 

i c 



(4) 



where /normal is the free-energy density in the normal state that is assumed to be independent of the 
magnetic field H = rot^, ^{f) is the complex order parameter, m* and e* are the effective mass and 
charge of the carriers [26j . We assume m* = 2m and e* ~ 2e, where m and e are the electron mass and 
charge, respectively, because the charge carriers in a superconductor are Cooper pairs. The coefficient 
a{T) changes sign at the superconducting transition T = Tc, while the coefficient f3{T) is always positive 
and depends on temperature only weakly (if at all). Equation ^ should be integrated over the volume 
of the superconducting sample (loop in our case) to obtain its free energy F. For simplicity, we will 
consider a circular ring with a diameter R and the same average perimeter as the square loops studied 
in our experiments. This will allow us to obtain at least a qualitative understanding of the experimental 
results without too complicated mathematics |27j . When the magnetic field H is perpendicular to the 
plane of the ring, i.e. directed along the z axis, A = e^Hp/2. In addition, if the thickness d of the ring 
is much smaller than the magnetic field penetration depth A(T), the magnetic field H can be considered 
equal to the external (applied) field. 

Minimizing the free energy F by variational methods leads to the famous GL differential equation 
for tp(f) ^6]. In Ref. [4] we have solved this equation numerically and have shown that this provides 
at least qualitative description of the experimental measurements. Here we will show that the problem 
can be treated analytically as well, provided that certain additional but reasonable assumptions are 
adopted. The central assumption is that both d and the width of the ring arm w are smaller than the 
GL coherence length £,{T). This allows us to neglect z- and p-dependences of the order parameter 
and write ip(f) = (</')• The free energy F is then minimized by ipni^ = iV'nl 6xp(m(/3) with integer n 
and [T5] 
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Here we assumed that w is much smaller than the average radius R of the ring and kept only the 
leading terms in w/R] cj) stands for ^/^o with $ = HttR^ the magnetic fiux through the ring and 
$0 = h/2e the magnetic fiux quantum. The state V'n(f) is called 'n-giant vortex state' because in this 
state, the ring hosts n magnetic vortices sharing the same core; we will also refer to it as 'n-vortex state' 
or simply 'state n' for brevity. The free energy Fn of the ring in the n-vortex state is: 
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where V = 2nRwd is the volume of the rmg and we introduced the GL coherence length f(T)^ = 
h^/Am\a{T)\. 

By noting that = ^f)Hc{T)'^ /2 with Hc{T) the thermodynamical critical magnetic field [SHj 

and by using the empirical expressions ^ H^{T) ^ Hc{0){l - t^), £_{T)^ = £_{0)^{1 + t^) /{I - t^) with 
t — T/Tc < 1, we can put Eq. ^ in the following form: 



where 7„ = [(,{0)/ R]'^[{n — (p)'^ + {ri^ /3 + 4>'^){w/2R)'^]. From the condition F„ = with ~ n we estimate 
the maximum number Umax of magnetic vortices that can be hosted by a ring to be nmax — ViR^ / w^{T) . 
For the rings with parameters corresponding to our experiments [Tc =1.18 K, ^(0) = 130 nm, w =230 
nm, R = 1.13 ^m), we find n^ax = 27, 37, and 43 at T = 1.00 K, 0.85 K, and 0.75 K, respectively. 

The heat capacity Cp of the ring in the n- vortex state is obtained by differentiating Fn twice with 
respect to the temperature: 



= = V^^^^^^^ [7^- 1 + 3^^(1 +7„n (8) 

We plot the heat capacity of the rings sXt — 0.95 (such a choice will be explained in the following) in 
different n- vortex states as a function of normalized magnetic flux in Fig.[S]by dashed lines [lines are the 
same in both panels] . When the magnetic field is varied continuously while the temperature is kept fixed, 
(as in our experiments), the ring can exhibit transitions between states with different n. In an increasing 
magnetic field, for example, the transition from n- to (n+ l)-vortex state becomes energetically favorable 
when Fn becomes larger than -Fn+i. The points 0„ where these (phase) transitions occur can be therefore 
foimd from the condition F.^ = -fn+i- This yields 0„ = (n + 1/2)(1 + u>^/12i?^), which appears to be 
exactly the points where Cp"^ and Cp"^^^ specific heat curves cross. In the thermodynamic equilibrium, 
therefore, the heat capacity Cp{(j)) should follow the lowest of the dashed curves in Fig. [51 thus exhibiting 
periodic variations with magnetic field with a period (pn+i — 4>n = 1 + /12R^ « 1, similar to what we 
observe in the experiments. 

A closer inspection of experimental results (Figs. [1] and |4] and Ref. [4]) reveals that in some of our 
measurements the peaks of Cp {H) appear more asymmetric than expected from the theory in the thermo- 
dynamic equilibrium, whereas in other measurements this asymmetry is not that striking. This suggests 
that at least in a part of our experiments the thermodynamic equilibrium is not reached. The rings 
exhibit a sort of non-equilibrium phase transitions, remaining in a state with a given number n of mag- 
netic vortices beyond the 'critical' magnetic field (/)„/7ri?^, where the thermodynamic transition to the 
{n + l)-vortex state should occur. To study such a possibility, we analyze the stability of the n-vortex 
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state ipnir) — |'0n| 6xp(m9j) with respect to an admixture of a fc-vortex state ipki^^ = IV'fel 6xp(iA;</?). 
Following Bezryadin et al. [7|, we write the free energy corresponding to the state ip^r) = V-'n(^) + '7'0fe(^) 
as 



F{t^) = Fn + Ajf + Bif (9) 

where A = V {iptf [a + 2/3 |V'„|^ + 7fc |a| ^(r) 7^(0)^] and B = (3\TPkf /2. In contrast to Ref. [7], we do 
not take the limit ^ in the above expressions and keep the terms of order [w/RY (but not the 
higher-order ones because our analysis is still limited to w <^ R). 

Because > 0, for A > Q the minimum of F{ri) is reached at = 0. This corresponds to the 
situation when the n- vortex state is stable and any admixture of a different, /c- vortex state is energetically 
unfavorable. In contrats, if ^ < 0, the minimum of Eq. ^ is reached for r; = ±(-^/2B)i/2 ^ 0. In 
this case the admixture of the fc-vortex state is energetically favorable and the n-vortex state becomes 
unstable. The condition A = Q defines therefore an instability threshold (or a 'superheating boundary' 
[3 [29]) of the state n. 

Let us now apply the above stability analysis to our experimental situation in which the magnetic 
field is slowly varied and transitions between different giant vortex states are observed. It is reasonable 
to assume that the transition from a state n to a state fc takes place when the former becomes unstable 
with respect to an admixture of the latter. Since in our experiments only transitions from n to n -I- 1 
(or to n — 1, in decreasing field) states are observed, we put k — n ± 1 and find that the rt- vortex state 
becomes unstable if (f> is increased above (decreased below) 



6± = nTl±\/2+[i?/e(T)]2 

/^rT^l2/9 2r./3T4/3 ]| ^^^^ 



y2RJ \ ^2+[R/aTW 



The instability of the n-vortex state in an increasing magnetic field occurs at larger fields than the 
thermodynamic transition does, i.e. (p^ > 0„. Similarly, in a decreasing field one has: i/)^ < 0„ (note 
that Eq. (|10p reduces to the result of Ref. [7] if we set u; = 0). This suggests that in an experiment 
the ring can remain in the state n (which in this case becomes metastable) beyond the thermodynamic 
transition and up to the instability point where a jump to n ± 1 state occurs. This is illustrated by solid 
lines in Fig. [6l These lines are obtained by following a given Cp"' curve up to (j)^ given by Eq. (fTO|) and 
then jumping to the next Cp"^^"* curve. The resulting dependence Cp{4>) is close to the experimental one 
(Figs. [Hand H]). On Fig. [6|the temperature was set tot — 0.95. A lower temperature would have lead to 
multiquanta transitions (see Ref. [4 ) . In our experiments - given a specific preparation of the sample as 
explained in the first section of this paper - we force the system to remain in the 10o oscillatory regime 
even at low temperatures. Thus to compare our signals with the model we used in Fig. [6|a temperature 
leading to a 10o regime. 
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Equation (|10p can also be used to study the periodicity of the Cp{H) dependence in more detail. 
The distance between consecutive instability points [jumps of Cp ((/>)] in the increasing magnetic field, 
A0+ = 0+ — 0n-i: sind in the decreasing field, Ai/),^ = 4'n+i ~ 'PVn ^'''^ found to be: 



We see that the periodicity of heat capacity jumps is not constant but changes with n. In particular, in the 
increasing field the distance between consecutive jumps decreases, whereas in the decreasing field, on the 
contrary, the distance between jumps increases. This conclusion can be compared with the experimental 
results shown in Figs. [3] and El For this purpose we plot the inverse of Ai/),^ as a function of n in Fig. 
[7] [30] . Since 1 / Ai/)^ roughly corresponds to the position of the maximum of the Fourier transform of 
Cp{H), the red solid line and the blue dashed line of Fig. [7]can be compared to Figs.[3]and[5]respectively. 
The predicted variation of the periodicity of heat capacity variations appears to be of the order of 10%, 
similarly to the experimentally observed values. Note that according to Eq. (fTT|) . the magnitude of this 
effect is controlled by the ratio w/R and that the effect disappears and the jumps of Cp{(j)) become equally 
spaced in the limit w 0. 

In the above analysis we put k = n ± 1, although nothing justifies such a choice. In general, solving 
A — for (p with arbitrary k yields a (rather lengthy) expression for the critical magnetic fluxes at which 
the 71-vortex state becomes unstable with respect to an admixture of the /c-vortex state. In the limit 
w/B. these are 



As the magnetic field modulus increases (decreases) we expect a transition to a state with k > n 
[k < n), and hence the '+' (' — ') sign should be chosen in Eq. (fT^ . The state k with respect to which the 
instability occurs first can be estimated by searching for a minimum (maximum) of 0,^^, as a function of 
k. This yields k — n± R/\/2^{T). We see therefore that at temperatures T close to Tc, when ^(T) oo 
and R/^{T) < 1, transitions to the states fc = n± 1 should indeed be favored. In contrast, at lower 
temperatures, when ^(T) becomes smaller than R, transitions to states k = n ±2, n ± 3, etc. can occur. 
Such 'multivortex jumps' have indeed been reported in Ref. ,4 and earlier in magnetic measurements of 
other groups [171 E] ■ They have also been analyzed theoretically [28] . This reminds us once again that 
these multivortex jumps or 'flux avalanches' are related to the existence of the metastable states and to 
the possibility for the system to remain in these states for a sufficiently long time. 

Another important comment is connected with the stability analysis presented above. As first noted by 
Horane et al. j29j . the stability analysis of a given n- vortex state ipnir) of a superconducting ring should 
not, in general, rely on a particular functional form [tpk{r) in our case] of the perturbation. By fixing 
the latter, we unavoidably overestimate the stability of the state under consideration. Unfortunately, the 
full treatment of the stability problem appears to be mathematically too involved and we therefore have 




(11) 



2n-k± ^2(n - kf + [i?/^(T)]2 



(12) 
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chosen to present a less general but much more physically transparent analysis here. Such an analysis 
still yield qualitatively correct conclusions. 

4 Discussion 

After the exposition of the experimental results, the data treatment and the development of the theoretical 
model describing our system, a deeper discussion of the main results is needed. Through the calculation 
of the SPWVD of the heat capacity signal versus magnetic field, it has been shown that the Cp{H) 
oscillations are not stationary. The periodicity of the jumps strongly depends on the sweeping direction 
of the magnetic field (increasing or decreasing). From Fig. 3 and Fig. 5, it is clear that when magnetic 
vortices are expelled from the loops (by decreasing the absolute value of the magnetic field), the heat 
capacity jumps are far from each other, giving a low- frequency oscillation of Cp{H). On the other hand, 
when the magnetic field is increased, vortices penetrate the loops and the jumps are closer and closer to 
each other, giving high-frequency heat capacity oscillations. This non-stationarity of the heat capacity 
oscillations is not negligible because from the experimental data, as well as from the theoretical estimation, 
the variation of the periodicity can be greater than 10% between expulsion of the first vortex and the 
penetration of the last one. 

In other words, and restraining our discussion only to positive magnetic fields, it can be said that as 
the field is increased, the system tends to spend less and less time in each successive giant vortex state. 
Furthermore when the vorticity n becomes high it appears from Fig. [6] that the occupied states are the 
stable ones, in contrast to the low field regime where the system almost always evolves along metastable 
states. This means that the higher the vorticity n, the easier it is to add another vortex into the loop, 
until the critical field is reached and the superconductivity suppressed. Thus the energy barriers the 
systems has to overcome to jump from the n to the n + \ state tend to disappear at higher fields. In 
that last case, the contour through which the fluxoid is quantized is close to the external edge of the 
loop and the supercurrents are localized near the outer boundary of the superconductor, leading to the 
" traditional" picture of the giant vortex state regarded as a surface superconductivity state [31] . 

On the other hand, in decreasing field, when the first vortices are expelled out of the loops (high 
magnetic field), the successive jumps occur with a large period, meaning that metastable states can 
survive for a long time before the energy barrier between the n and n—1 states is suppressed. This feature 
is confirmed by Fig. [6l near the critical field the system spends more time in metastable states when 
vortices are expelled than it does in increasing field. In this situation the apparent contour through which 
the fluxoid is quantized is close to the inner boundary of the loop; the kinetic energy of the supercurrent 
carriers is higher than it would be if the system was in its ground state, where the supercurrents are 
localized at the outer boundary. 

Furthermore the differences between penetration and expulsion of vortices are enhanced due to the 
fact that in increasing field the surface barrier can be destroyed by surface defects, which is not the case 
in decreasing field [32j • This effect has not been taken into account in the theoritical work leading to Eq. 
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but it is in good agreement with the smaU deviation of the experimental data presented in Fig. [3] 
from the hnear behavior predicted by Eq. [TT] and Fig. [71 

5 Conclusions 

We have studied the heat capacity behavior of an assembly of non-interacting superconducting loops in 
applied magnetic field. Heat capacity discontinuities are observed when the magnetic field is swept. The 
jumps are associated to vortex expulsion or vortex entrance in the loops. We have studied the variation 
of the periodicity of the heat capacity jumps versus the magnetic field. Using the signal processing tool 
based on Wigner-Ville distribution, we were able to image the variation of that periodicity versus the 
direction of the magnetic field sweep. As the field is increased, the vortices entering the loops are not 
submitted to the same barrier of energy as when the field is decreased and the vortices are expelled. The 
GL theory describes our experimental observations with reasonable accuracy. 
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H (mT) 

Figure 1: Heat capacity measurement of the superconducting rings at 0.75 K. The heat capacity jumps 
at -35 mT and 35 mT corresponds to the transition between the superconducting state and the normal 
state. The insets focus on two areas after signal processing reducing the visible noise (10 points adjacent 
averaging): small heat capacity oscillations appears to be superimposed on a slowly varying background. 
The arrow indicates that in this experiment the field is swept from -45 mT to 45 mT. 
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Figure 2: Modulus of the Fourier transform of the signal shown in Fig. 1 windowed from H=-30 mT to 
H=30 mT (superconducting part of the plot), after low pass filtering to remove the background. A large 
peak appears at the frequency v = 1.65 mT~^, which corresponds to one magnetic flux quantum in a 
square of 1.86 /um side. The inset shows an electron micrograph of one loop : the external side is 2.05 
/um and the internal side is 1.65 /xm. 
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Figure 3: (Color online) This smoothed pseudo Wigner-Ville distribution (SPWVD) is calculated from 
the signal presented in Fig. 1 after low pass filtering to remove the slowly varying background. The 
upper part of the graph shows the heat capacity signal after substraction of that trend. The lower part 
is the SPWVD of the signal (a magnetic field- frequency representation of the heat capacity data). The 
evolution of the periodicity versus the magnetic field appears clearly in the middle of the graph, where 
the frequency v of the oscillations varies from 1.50 mT^^ to 1.78 mT^^, which represents more than 
a 15% change. The color code of this graph is only used to illustrate the difference in the calculated 
amplitude: from blue for the lowest amplitudes to red for the highest ones. 
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Figure 4: Heat capacity measurement of superconducting rings at 0.75 K in decreasing magnetic field. 
The oscillations with a periodicity of one (po are expanded in the inset. The arrow indicates that in this 
experiment the field is swept from +45 mT to mT. 
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Figure 5: (Color online) This SPWVD is calculated from the signal presented in Fig. [4] after low pass 
filtering to remove the slowly varying background (upper part of the graph). The lower part is the 
SPWVD of that signal. The evolution of frequency versus magnetic field is inverted as compared to the 
case where the magnetic field is increased (see Fig. [3]). The color code is the same as in Fig. [31 
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Magnetic flux {|) = O/Oq 

Figure 6: (Color online) Heat capacity of a superconducting ring calculated from the Ginzburg-Landau 
theory at T/Tc = 0.95. The ring is assumed to be circular with the same average perimeter as the square 
loop studied in the experiments and with the same arm width and thickness. Dashed lines correspond 
to giant vortex states with vorticity n. Transitions between consecutive n-vortex states taking place at 
'critical' magnetic fluxes given by Eq. (jlOp . yield the (red) solid curve in the upper panel in the increasing 
magnetic field and the (blue) solid curve in the lower panel in the decreasing magnetic field. The large 
jumps of the heat capacity at ^ ~ 14.5 correspond to the normal-superconducting transition. 
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Figure 7: (Color online) Inverse periodicity of the heat capacity Cp{(f>) [Eq. (|lip ] at T/Tc = 0.95 as a 
function of the magnetic vortex number n. The inverse periodicity increases with n in the increasing 
magnetic field (red solid line) and decreases with n in the decreasing magnetic field (blue dashed line). 
The tendency is the same as found experimentally (Figs. [2]and[5l respectively). The magnitude of the 
effect ^ 10% is also comparable to the measured one. 
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